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ABSTRACT 

Recent coronagraphic imaging of the AB Aurigae disk has revealed a region of low polarized scattered 
light suggestive of perturbations from a planet at a radius of ^ 100 AU. We model this darkened 
region using our fully non-plane-parallel radiative-transfer code combined with a simple hydrostatic 
equilibirum approximation to self-consistently solve for the structure of the disk surface as seen in 
scattered light. By comparing the observations to our models, we find that the observations are 
consistent with the absence of a planet, with an upper limit of 1 Jupiter mass. 

Subject headings: accretion, accretion disks — circumstellar matter — stars: individual (AB Aur) — 
planetary systems: protoplanetary disks — planetary systems: formation 
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1. INTRODUCTION 

The Herbig Ae star, AB Aurigae (AB Aur), hosts 
a much-studied disk of gas and dust thought to be 
representative of protoplanetary disks during the Jo- 
vian planet formation phase. The circumstellar ma- 
terial around AB Aur has been imaged in scattered 
l ight by several ob servatorie s, including STIS on HST 
(jGradv et al .l llQQQl). Subaru (iFukagawa et al.l l2004D. and 



the Lyot Project ( Oppenheimer et al.l |2008[ ) . These ob- 
servations show that AB Aur is surrounded by what ap- 
pears to be a nearly face-on disk, with structures that 
look like spiral arms. 

Using adaptive-optics c oronagraphy and polarimetry, 
lOppenheimer et all ()2008f ) imaged the AB Aur disk be- 
tween radii of 43 — 302 AU. Their images revealed an 
azimuthal gap in polarized light at a radius of about 100 
AU, along with a 2-a bright "spot." They interpreted 
this to indicate the presence of a possible massive planet 
of 5 — 37 Jupiter masses (Mj). 

The position angle (PA) of the putative planet is coin- 
cident with the minor axis of the disk, as a well as a gap 
betwee n two spiral arms obser ved in other scattered light 
images (Fukagawa et al. 2004), suggesting that the dark- 
ened region results from anisotropics in the overall disk 
structure, rather than locally confined to a single per- 
turber in the d isk. On the other hand, t he spiral arms 
are not seen by lOppenheimer et"al] ()2008l ) . 

Another alternative is that there is no true gap in the 
disk at all, but rather that the lack of scattered polar- 
ized light in this region is purely a geometrical effect 
caused by the inclination of the disk (jPerrin et al.| [2009). 
In this scenario, the disk is inclined enough that the far 
edge of the disk is back scattering rather than forward 
sacttering, creating a region of lower total polarized in- 
tensity in the back scattering region. This is verified by 
a lower polarization fraction in the supposed gap, but 
no corresponding decrease in the total scattered light. 
This is likely to be the correct interpretation of the AB 
Aur observation. The objective of this paper is to put 
firm limits on the mass of any planet embedded in AB 
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Aur, complemetary to the interpretation bv lPerrin et al.l 
(2009). 

Planets embedded in optically thick accretion disks, 
like the disk around AB Aurigae, are expected to pro- 
duce perturbations in the density and temperature s truc- 
ture of the disk. iJang-Con dell fc Sasselovl ()2QQ3l ) and 
iJang-Condell Sasselovl (12004) calculated the magni- 
tudes of these perturbations for a range of planet masses 
and distances. They predicted the formation of a shadow 
at the position of the plan et paired with a brig htening 
just beyond the shadow. iJang-Condelll (|2008L hence- 
forth Paper I) improved upon these calculations, by 
self-consistently calculating the temperature and density 
structures under the assumption of hydrostatic equilib- 
rium and taking the full three-dimensional shape of the 
disk into account rather than assuming a plane-parallel 

disk. 

lOppenheimer et"al] (|2C 



suggested that the observed 
structure resembled models of dust t rapped in mean mo- 
tion resonances with a planet (e.g. 'Kuchner fc Hol manI 
l2003f ). However, AB Aur hosts a gas-dominated disk, as 
confirmed by detectio ns of numerous molecular s pecies 
(iMa nnings fc SargentI I1997L 120001: iDent et all 120051: 
Pietu et all 120051 : iCorder et al.l 120051 : iSchrever et al.l 
2008) In such a disk, the dynamics of the gas domi- 
nate the system and the dust traces the gas on orbital 
timescales. Moreover, the disk is optically thick, so the 
scattered light image traces only features in the upper 
surface of the dis k. Thus, the darkene d region and bright 
"spot" imaged by Oppenh eimer et al.l (|200 8) are unlikely 
to trace dust concentrations, as might be the case were 
the disk optically thin, but could perhaps be scattering 
off structure in the disk created in the wake of a planet. 

Here we model the observed structure as a shadow 
caused by the tug of a planet in approximate hydro- 
static equilibrium with the gas disk, rather than reso- 
nant trapping of ballistic grains. This interpretation is 
consistent with short stop ping time of the dust g rains 
in the gas, - 10"^ orbits (jYoudin fc Chiang|l2004f ). We 
use the algor i thms and code described in Paper I and 
iJang-Condelll (|2009L henceforth Paper II) to model the 
AB Aurigae system (stellar mass 2.4±0.2 M©, luminosity 
48 1/( 7), accretion rate 10~^ M0yr~^: iGarcia Lopez et all 
120061 : iRodgersI l2001f ) . We synthesize images of the disk 
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in scattered light and compare them with the corona- 
graphic images of the AB Aurigae disk to constrain the 
mass of the perturbing planet. 

2. MODEL DESCRIPTION 

The model we adopt is described in detail in Papers I 
and II, with parameters adjusted for the AB Aur system. 
The stellar parameters we adopt are mass = 2.4 M©, 
radius = 2ARq, and effective temperature Teff = 10^ 
K, which are consistent with L* = 481/0. For the disk 
model, we assume a viscosity parameter of a = 0.01 and 
accretion rate M = lO~^M0yr~^. Modifying these last 
two parameters effectively changes only the overall mass 
of the disk. Since we are only examining the perturbative 
effects of hypothetical planets in the disk, changing a or 
M do not significantly alter our final conclusions. Indeed, 
adopting M = 10~^ M0yr~^ was not found to alter our 
overall results. 

The amount of heating from stellar irradiation in the 
disk is sensitive to the opacities used to calculate ra- 
diative transfer in the disk. Dust opacities are perhaps 
the least observationally constrained parameter for disk 
modeling. Hence, we calculate the opacities usin g a Mi e 
scattering code developed by iPollack fc Cuzzil (119801). 
The c omposition of the dust is that used in IPollack et al.l 
(|1994f ). a mixture of water, troilite, astronomical sili- 
cates, and organics. We adopt a size distribution for 
the dust of of n{a) oc a~^-^ where a is the radius of 
the grain, with maximum/minimum radii of 1 mm/0.005 
microns. The large maximum grain size represents co- 
agulation of grains in protoplanetary disks. The dust 
is well-mixed with the gas and uniform in mass opacity 
throughout the disk. The mean absorption and total ex- 
tinction (absorption+scattering) averaged over a black- 
body of 10,000 K representing the stellar spectrum are 
tvp = 2.68 and Xp = 17.3 cm^-g~^, respectively. For ra- 
diation at disk temperatures, the relevant opacities are 
the flux-averaged absorption and Rosseland mean opac- 
ity, i^R = 2.73 and hzp = 1.46 cm^-g~^, respectively, 
which are calculated for a blackbody at 100 K. The tem- 
perature at the surface of the disk in the resulting models 
ranges from 64 to 143 K over the computation domain, 
so a choice of 100 K is a reasonable compromise. 

Scattered light dominates the surface brightness of the 
disk at 1.64 /im. We calculate scattered light as described 
in Paper II. The idealized images have been convolved 
by a Gaussian PSF with FWHM 0.094'', the diffraction- 
limited resolution of the AEOS teles cope, to compare to 
the Lyot Project images of AB Aur (jOppenheimer et al.l 
|2008). 

To compare our model to the "Oppe nheimer et all 
(l2008l) P image in scattered polarized light, we must in- 
clude the effects of polarization. The size distribution of 
the dust grains is of n{a) ex a~^"^, so scattering is domi- 
nated by the smallest grains. Thus, Rayleigh scattering 
is a reasonable approximation for the angle-dependent 
polariz ation caused by scattering from dust grains (e.g. 
iGraham et al.ll2007l ). The f ractional linear polariz ation 
from Rayleigh scattering is (|Hansen fc TravisllTOTl 
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where is the angle of deflection. We assume that cir- 
cular polarization is negligible, so V = 0. Here, /, Q, 
and V are the canonical Stokes parameters. 

In our model, the scatterers are dust grains in the layer 
of the disk that is optically thin to stellar light. We 
assume isotropic single scattering in this optically thin 
layer of the disk to model the observations. Although 
our model incorporates thermalization caused by multi- 
ple scattering of stellar photons, it does not include the 
effect of multiple scattering on the polarization of the 
scattered flux. Monte Carlo calculations of multiple scat- 
tering in circumstellar disks (Wood et al.l ll996f ) suggest 
that multiple scattering can increase the net polarization 
from a smooth disk by a factor of 1.5 to 2. But the issue 
that is important to our study is not the net polariza- 
tion, it is whether the dimple caused by a planet would 
be enhanced by multiple scattering effects; we will leave 
this second-order effect to a future study. 

3. COMPARING THE MODELS TO THE DATA 

The simulated images of the model disk depend on 
the assumed inclination and orientation. This has been 
variously reported for AB Aur as 12 — 30° with PA 
60-80° ([Marinas et al.l l2QQ6D. 25 - 35° with PA 50 - 60° 
( Fukagawa et all I2004D . 23 - 43° with PA 58 - 63° 
(iPietuet all 120051 ). and 21° with PA 59° (jCorder et all 
'2005'). Suffice it to say that the inclination of the disk 
is highly uncertain, although the PA is largely agreed 
to be around 60°. We adopt an inclination of 25° as 
being a reasonable compromise between the various ob- 
servations. For simplicity, we assume that the PA of the 
supposed planet is equal to the PA of the minor axis of 
the inclined disk. 

Fig. [1] shows the real and simulated P images of AB 
Aur in scattered polarize d light at 1.64 micr ons, with 
the observed image from I Oppenheimer et all ([2008) in 
the lower left. We examine disk models with planets 
of 0, 100, 300, and 900 (0, 0.33, 1, and 3 Jupiter 
masses), as indicated. The color scale has been calibrated 
to show the same dynamic range in both the data and the 
model images. The blacked out inner circle represents the 
coronograph in the observed image, and the simulation 
boundary in the model images. The northwest edge of 
the disk tilted away from the observer. 

There is some amount of uncertainty as to the exact 
location of the star on the observed image, while in the 
simulat ed models the stellar posi tion is known exactly. 
Indeed. lOppenheimer et al.l (|2008l ) report an offset in the 
photocenter of the inner disk of about 88 mas or 12 AU 
along the minor axis. In order to make a fair comparison 
between the data and model, we recenter the images by 
the following procedure. 

We select an isophotal contour, with the brightness 
level set by the bright "spot" in the data image. Al- 
though we will henceforth refer to this position as the 
"planet," we understand it to represent an image feature 
that may or not correspond to an actual planet. For 
the model images, we determine the local maximum in 
brightness in the dimple caused by the 1 Jupiter mass 
planet, and use this brightness level across all model 
images. The isophotes are indicated in Fig. [T] by solid 
black lines. We then calculate the best fit ellipse to this 
isophote, indicated by the long-dashed white lines. The 
remaining calculations are now all calculated using the 
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Fig. 1. — Observed P image and simulated images of AB Aur, scaled to the same spatial resolution and dynamic range. The observed 
image is in the lower left. The blacked out inner circle represents the coronograph in the observed image, and the simulation boundary 
in the model images. The disk models contain planets of (top center), 100 (top right), 300 (bottom center), and 900 (bottom right) 
M0. The long-dashed while lines show the best-fitting ellipse to the isophotes indicated by black lines. The solid lines are drawn from the 
center of the ellipse through the position of the planet, which is equal to the minor axis in the model images. The minor axis is indicated 
by a black line in the observed image. The dotted white lines show the minor axis opposite the planet. The dashed and dot-dashed lines 
indicate the major axes of the ellipse. 



center of the ellipse as the image center. The calculated 
PA for the data is 55°, while the planet's PA is -25°. 
The angle between minor axis and planet PA is small, 
so aligning the planet with the inclination vector of the 
model disks is a reasonable assumption. 

The major and minor axes of the best fitting ellipses 
for the observed and simulated images are tabulated in 
TablelH assuming that V = 144 AU. The offset listed for 
the model images is the distance between the star's posi- 
tion and the center of the ellipse. The best-fit ellipses for 
the model images are slightly more eccentric than that 
for the data, indicating that either the inclination of 25° 
is too high, or that the scattering properties are incom- 
pletely modeled. For example, a small amount of forward 
scattering or slight flattening of the flared disk structure 
could both create a less eccentric best fit ellipse. How- 
ever, we are interested in studying the perturbative effect 
of a planet on the overall disk structure, not creating an 
exact match to all the disk properties. Indeed, changing 
parameters such as disk inclination, accretion rate, and 
disk mass do not significantly change our overall results 
on the upper limit for the mass of a planet. 

In Fig. [TJ the white solid, dashed, dotted, and dash- 
dotted lines mark the line through the planet, the north- 
east major axis, the southeast minor axis, and southwest 
major axis, respectively. In the simulated images, the 
planet is aligned with the minor axis, and only one ma- 
jor axis cut is displayed because of symmetry. In the 
observed image, the minor axis is marked by a black line 
in the lower left of Fig. [H offset from the solid white line 



TABLE 1 

Best-fit ellipse parameters to isophotal contour, in AU. 





major axis 


minor axis 


offset 


xi 


data 


129 


108 






model: 










no planet 


131 


104 


31 


0.76 


100 Me 


131 


103 


31 


0.44 


300 Me 


132 


100 


29 


2.18 


900 Me 


133 


97 


26 


16.8 



through the planet. 

In Fig. [21 we plot the surface brightness profiles along 
the cuts indicated in for the observed and simulated im- 
ages, with the solid, dotted, dashed, and dot-dashed lines 
indicating the profiles along the similarly marked cuts in 
in Fig. [TJ with the observed profiles in the lower left, 
and the remaining plots corresponding to the indicated 
planet masses. In the observed data, we plot an addi- 
tional grey line, the profile along the northwest minor 
axis, indicated by the black line in Fig. [H This is to 
show that the profile does not deviate significantly from 
the profile along the plane t line, except for the bright 
"spot" pointed out in (Op penheimer et al.ll2QQ8T ). The 
distances of the major and minor axes of the elliptical fit 
to the isophotes are indicated by down and up arrows, 
respectively. 

We calculate the error bars according to the following 
procedure. At a given radial distance from the ellipse 
center, we sample the brightness along the circle of the 
given radius. The inclination of the disk introduces an 
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Fig. 2. — Surface brightness profiles along different cuts through the observed and simulated disk images. The observed data are plotted 
in the lower left, while the remaining plots show the profiles for models with planets of (top center), 100 (top right), 300 (bottom center), 
and 900 (bottom right) M0. The lengths of the major and minor axes if the isophotal ellipse are indicated by down and up arrows, 
respectively. The solid lines are the brightness profiles from the center of the ellipse through the position of the planet. The grey line in 
the lower left plot is the profile along the minor axis of the isophotal ellipse. The dotted lines are profiles along the minor axis opposite 
the planet. The dashed and dot-dashed lines are profiles along the major axes of the ellipse. 



intrinsic variation in brightness that is symmetric about 
only on axis, so we find the best-fitting function of the 
form a + bsm{0 + c) + dsm{20 + e) where is the po- 
sition angle and calculate the rms deviation from this 
best-fitting function to be the error. The errors are sim- 
ilar in magnitude to those calculated with a best-fitting 
function of a -\- b sm{6 + c) so we are not over-fitting the 
data. 

The surface brightness profiles of the data (Fig. [2] lower 
left) are qualitatively simular to those of the planet-less 
model (Fig. [2] upper middle), with the steepest profile 
along the southern minor axis (dotted line), and the 
northern minor axis (solid line) consistently dimmer than 
the major axes (dashed and dot-dashed). This validates 
our assumption that the northwestern edge of the disk is 
tilted away from the observer. 

When a planet is added to the image, it creates a dim- 
ming in the surface brightness profile slightly inward of 
the length of the minor axis and a brightening just out- 
ward of it. The magnitude of this S-shaped perturbation 
grows with increasing planet mass. The shifts in the pro- 
files along the other axes can be explained by the shift 
of the ellipse center as the planet mass increases. 

4. UPPER MASS LIMIT 

In order to establish how well the models match the 
data, we need to quantify the goodness of the fits. Since 
we are interested in the perturbations caused by the pres- 
ence of a planet and only qualitatively interested in the 
rest of the disk structure, we re-scale the brightness pro- 
files before directly comparing the data to the models. 



For a given disk image, we first scale the radius to the 
length of the minor or major axis, as appropriate, to 
derive a normalized distance. In the data image, the 
elliptical radius along the planet's position differs from 
the minor axis by 1%, so it makes little difference to use 
the minor axis to scale this profile. We then divide the 
brightness profile along the planet axis by that of the 
major axis. This cancels out any small variations that 
might result from uncertainties in inclination angle. In 
the case of the data image, we average between the pro- 
files of both major axes before scaling. 

The resulting normalized brightness profiles for the 
data and models are plotted in Fig. [3l Using these pro- 
files and error bars measured from the data, we calculate 
the reduced over the normalized radius range from 0.6 
to 1.3. The resulting values of reduced xi tabulated 
in Table [TJ We find that both a Saturn mass planet and 
no planet at all are good fits to the data, while a Jupiter 
mass or larger is excluded. 

5. CONCLUSIONS 

We find that scattered polarized light images of AB 
Aur do not indicate the presence of a massive planet. We 
put an upper limit of 1 Mj on any planet at this location, 
well below the s u^^ested mass of 5 — 37 Mj postulated by 
i Qppenheimer et al. (2008). Indeed, the best models are 
consistent with t he absence of any planet, and support 
the hypothesis of iPerrin et al.l (|2009l ). that the "gap" is 
simply an inclination effect rather than an indication of 
a planet. 

On the other hand, if a massive planet several times 
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Fig. 3. — Scaled surface brightness profiles through the planet 
position. The data are plotted as diamonds with error bars. Pro- 
files of disk models with 0, 100, 300, 900 planets are shown 
with solid, dotted, dashed, and dot-dashed lines, respectively. 



Jupiter's mass did exist in the disk, our models suggest 
that it would have been detectable. In other words, our 
work suggests that massive planets on 100 AU orbits can 
now be detected via scattered light imaging of nearby 
Herbig Ae/Be disks. 

As a Herbig Ae star, AB Aur will eventually evolve 
into an A- type main sequence star. We can compare AB 
Aur to two main sequence A-type stars with directly im- 
aged planets: Fomalhaut ' s plane t is at 120 AU separation 
from the star (|Kalas et al.ll2008[ ). an d HR 8799 's planets 
are at 24, 38, and 68 AU separation (iMarois et al. 2008). 
Assuming no significant planetary migration, our upper 



limits indicate that AB Aur is not a younger analog of 
Fomalhaut. The coronographic spot covers the inner 40 
AU of AB Aur's disk, so it may yet prove to be an ana- 
log of HR 8799. If planets do migrate significantly after 
dissipation of the gas disk, and there is some mechanism 
for expelling planets out to large distances, then AB Aur 
may evolve into a Formalhaut analog. 

If systems like Fomalhaut and HR 8799 are common 
and planet formation and migration occur early, then 
there may be many planets waiting to be detected in 
known Herbig Ae/Be disks. The requirements include 
high angular resolution (< 0.1^'), good star light suppres- 
sion, and small inner working angle (< 0.3^'). Ongoing 
coronographic surveys, like the SEEDS survey on Sub- 
aru and future surveys on the Keck, VLT, and Gemini 
telescopes, will have the power to probe this regime. 

The authors are indebted to Ben Oppenheimer for gen- 
erously providing the observational data needed for this 
paper. We also thank an anonymous referee for helpful 
comments that greatly improved our paper. Simulations 
presented in this paper were carried out using the "borg" 
cluster administered by the Center for Theory and Com- 
putation of the Department of Astronomy at the Univer- 
sity of Maryland. Resources supporting this work were 
also provided by the NASA High-End Computing (HEC) 
Program through the NASA Center for Computational 
Sciences (NCCS) at Goddard Space Flight Center. HJ- 
C acknowledges support from a Michelson Postdoctoral 
Fellowship under contract with the Jet Propulsion Lab- 
oratory (JPL). JPL is managed for NASA by the Cali- 
fornia Institute of Technology. 
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